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1 Introduction 

Deriving the equations of motion of a test particle from field equations of general relativity has 
a long history, see, for example, the account in [1]. It is now well established in the literature 
that the equations of motion need not to be postulated separately, but can be derived from 
the field equations of general relativity. More precisely, the conservation of the ordinary 
matter energy-momentum tensor, i.e. V^T'^'^ = 0, is enough to find the geodesic equation 
and also the equations of motion of a spinning test particle (pole-dipole particle). In general 
relativity, a pole-dipole particle is described by the four-momentum and the tensor of spin, 
and the dynamics is governed by the Mathisson-Papapetrou equations [2, 3]. Combining these 
equations with a proper supplementary condition, one gets a self-consistent set of equations 
for describing the motion of a pole-dipole particle in a given space-time, for example see 
[4-8]. The case in which both gravitational and electromagnetic fields are present and the 
particle is charged was studied by Dixon and Souriau [9, 10]. The equations of motion of a 
charged spinning test particle are known as Dixon-Souriau equations. For some solutions of 
the Dixon-Souriau equations in the specific backgrounds, see [11-13]. 

The multipole approximation method of [2, 3] has been used to investigate the motion 
of a test particle in the context of some modified gravity theories, for example see [14-16] . In 
this paper, using the multipole method, we derive the equations of motion of a pole-dipole 
test particle in modified gravity theories in which the covariant divergence of the ordinary 
matter energy- momentum tensor is given by V^T'^'^ = A'^ , where the nonminimal coupling 
term A'^ is an arbitrary covariant vector, and in each modified gravity theory this vector 
can be determined exactly when the particular way by which the gravitational fields of the 
theory couple with the ordinary matter is specified. It is obvious that A'^ leads to a modi- 
fication of the equation of motion of the test particle. In fact, we generalize the Mathisson- 
Papapetrou equations in such a way that they can be used in any torsionless non-metric 
theory of gravity. It should be noted that, in metric theories of gravity, the ordinary matter 
energy-momentum tensor is conserved, and consequently the single-pole particles motion is 
described by the geodesic equation and the pole-dipole particles motion is governed by the 
Mathisson-Papapetrou equations. Furthermore, we make use of the generalized Mathisson- 
Papapetrou equations in order to investigate the test particle equations of motion in Modified 
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Gravity (MOG) [20]. In addition, we derive some conserved quantities for the test particle 
motion. 

2 Energy-momentum conservation 

Let us study the test particle motion in a modified gravity theory without torsion in which 
the usual conservation law of the ordinary matter energy-momentum tensor is replaced by 

V^T^^ = (2.1) 

We know that modified gravity theories which yield to such an identity for the energy- 
momentum tensor, would violate the Einstein equivalence principle. In fact, one can show 
that the covariant divergence of the energy-momentum tensor is zero if the gravitational 
fields (such as scalar fields or vector fields) of the theory do not interact directly with the 
matter. In the presence of the test particle, the components of the metric tensor are g'^A = 
9fj,x + ^QplXi where g^j^x is the background metric and 6g^\ is the metric perturbation caused 
by the test particle, and it is supposed to be small. Also, let the energy-momentum tensor be 
ji^A _ J<^lX _|_ ^j'M-*' and the nonminimal coupling term be A'^ = A'^ + 6A^, where the tensor 
jjT/iA (jggcj-ibgg the material distribution inside the test particle and 5A'^ stands for the test 
particle's contribution to the coupling term. If we assume that the test particle's trajectory 
is outside the massive bodies then the background energy-momentum tensor T^^ vanishes 
inside as well as near the test particle. Thus, using the identity (2.1) we get 

V^r/^^ = 5f^ (2.2) 

where t^'^ = \J—g ST^^ and Sf"^ = \/—g SA"^. As we mentioned before, the use of (2.2) 
alone is sufficient to derive the equations of motion of the test particle, and we do not 
need the filed equations of the theory. In the following sections by using this equation we 
derive the generalized version of the geodesic equation. Furthermore, we find the generalized 
Mathisson-Papapetrou equations for describing a spinning test particle. 

3 The single-pole particle 

In this section, we derive the equations of motion for a single-pole test body. Our analysis is 
based on the integration of the conservation law (2.2) over the world tube of the test body. 
This procedure is independent of a specific choice of energy-momentum tensor for the test 
particle. In fact, in the four-dimensional diagram the interior of the particle can be considered 
as a tube-like region (world tube) contained in the 3-dimensional time-like hypersurface S(i). 
A representative continuous curve through the tube is parametrized by X^{t). Coordinates 
within the word tube with respect to a coordinate system centered on X'^(t) are labeled by 
x^, see [3] for details. 

Equation (2.2) can be rewritten as 

d,T>^Hx) + r;;„(x)T'^"(x) = 6f\x) (3.1) 

where r^„(x) are the Christoffel symbols corresponding to the background metric Since 
gf^i, changes very little inside the test particle, the Christoffel symbols can be developed in 
Taylor series around the point X'^{t) 

^Li^) = ^UX) + 5x^d^TUx) + ... (3.2) 
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where Sx'^ = x"' — X'^. The single-pole particle has the simplest internal structure and its 
dipole as well as higher multipole moments are zero 

[ dx^T^'^dPx = 0, f dx^dx'^T^^d^x = 0, ... (3.3) 

Note that the integrals are taken over the 3-dimensional hyperspace which denotes the 
interior of the test particle at time t. Substituting (3.2) into (3.1) and taking into account 
that we are considering a single-pole particle, we get 

d,T^\x) + r;;jX)T^"(x) = 5f\x) (3.4) 

Integrating this equation over the hypersurface S and taking into account that t^'^ are zero 
outside S, we find 



— j T^'^d^x + r;„(X) j r^"(i^x = j 6rd-'x (3.5) 
On the other hand, multiplication of equation (3.4) by x'^ yields 

t''^ = d^ix^T"^) + x''T^^^{X)t^" - x^df^ (3.6) 
Integrating this equation over S, we derive 

J T^^cfx = j^j x^T^^Sx + Vl^{X) j x^r^^d^^ - j x''df^(fx (3.7) 

Furthermore, using the requirements (3.3), we can write 

J x^rf'^cfx = X^it) J Tf'^cfx (3.8) 
j x^df^d^x = X^it) j 6f^d^x 
Substituting these equations into (3.7) and using equation (3.5), one gets 

r-^d'x = ^ / r^°d3x (3.9) 



dt 

With the help of the A = component, we rewrite (3.9) as follows 



dt dt 

Combining equations (3.5) and (3.10) we find 

d fdX^ f nn \ ,^^^dXf'dX' 



r-^d\='^'^ I r'^d'x (3.10) 



r'^d^x + rUx)--—- / rOOd^x = / 6f'd^x (3.11) 



dt\ dt J / ^"'^ ' dt dt 



On the other hand, the line-element of the curve X'^lt) is given by ds'^ = g^^ydX^dX'^ . Also, 
the corresponding four-velocity is = Thus equation (3.11) becomes 

^(mu^) + mr^^n^u" = j 6f^d^x (3.12) 
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where m is defined as follows 

1 ' .00 j3 



m= y T^Vx (3.13) 
Multiplying equation (3.12) by ux and taking into account that u'^'V^ua = 0, we find 

^ = ^0^;^ f 5fd^x (3.14) 
as J 

In general relativity, the right hand side of this equation is zero and so m is a constant of the 
motion. In this case m is interpreted as the rest-mass of the particle [3]. Substituting (3.14) 
into (3.12) we find 

here we introduced A°" as 

A'' = u° I dfd^x (3.16) 



m ( ^ + r^„n^n° ] = {6^ - ^x^u^)A'^ (3.15) 



Finally, equation (3.15) can be brought into a more compact form, namely 

= {6^ - u\^)A'' (3.17) 

where ^ = u^V ^ is the directional covariant derivative. Equation (3.17) is the generalization 
of the geodesic equation. As we expect, orbits of single-pole particles are not the geodesies of 
the background metric. Prom (3.17), one can immediately read off the additional contribution 
to the equations of motion due to the nonzero covariant divergence of the matter energy- 
momentum tensor. One can interpret the left-hand side of (3.17) as an extra force - "extra" 
in comparison to the case A'^ = - influencing the test particle motion. Such an extra 
force due to the nonminimal coupling term on the right-hand side of (3.17) may provide an 
interpretation for the observed mass discrepancy in spiral galaxies and clusters of galaxies 
[17, 18], or in the context of the so-called Pioneer anomaly [19]. Por some attempts to 
construct consistent modified gravity theories that might yield to an appropriate extra force 
for explaining the aforementioned problems, see [20, 21] and references therein. 



4 The pole-dipole particle 

In this section, we derive the equations of motion of a pole-dipole particle. The pole-dipole 
particle is a particle for which all integrals J dx'^^Sx^'^ ...5x^"T^'^d^x with more than one 
factor 5xP vanish. Remembering that the test particle is a pole-dipole particle, we substitute 
equation (3.2) into (3.1) and then integrate over the hypersurface S. The result is 



d 



j T^^d^x + T^^^{X) j T'"^d^x + d^Tf^^{X) J fe''r^"d3x = J 5f^d^x (4.1) 
Also, multiplying equation (3.1) by x'^ and integrating over S, it follows that 

j r^0d3x + r^„(X) j dx^T'^'d^x- j Sx^df^d^x (4.2) 



'r'^^d'x = ^ [ 6x^T^°d^x- '^^'^ 



dt / dt 
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Here we have used (4.1) to simplify (4.2). Another useful equation can be obtained by 
multiplying equation (3.4) by x°'x^ and integrating over S. The result is 

|,..,.«.3, , ^ ^ / . / A (4.3) 

Equations (4.1)-(4.3) contain enough information to describe the motion of pole-dipole par- 
ticles. In the following, we will bring these equations into more familiar form. To do so, 
analogously to [3] we introduce the quantities M^"^, S""^ and M"^ as follows 

M^"^ = -n° J dx^T^^cfx (4.4) 
M"^ = n° / T^'^d^x (4.5) 



= - M/^°°) (4.6) 

where S"'^ is the total angular momentum or macroscopic spin of the test particle. Further- 
more, one should note that M^°^^ = (because 6x^ is zero) and M""'^ = —vPS°^^. Also, we 
define the new quantity A"'^ as follows 

= u^ j 5x'^5f^d^x (4.7) 
With these definitions the equations (4.1)-(4.3) become 



ds 



( ^ 1 + r;:,M^° - 9,r;:„M'^^" = (4.8) 



M'^^ = — - 3^ ( —IT - r^a^''''" - (4-9) 



Ifs 

0(^ja/3A ^ j^PaX^ ^ u^^M^^^ + n^M"^° (4.10) 



Cyclic permutation of the indices in equation (4.10) and subtraction of the second from the 
combination of the first and the third of the permutations yields 



By setting A = in equation (4.9), we find 



M'^o = ^M°° - - r°^M^^" - A'^o (4.12) 

If we contract (4.12) with and use the same choice for the mass as in [3], namely 

1 

we obtain 

m = ^ (mOO + rO,5^%'^) + - ^A'^o (4.14) 



V 



m= — (M'^o + r!;„n^5'^°) (4.13) 



(n0)2 ^ ' ' 'vP ds 
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This mass parameter sometimes called the " kinematical" or "monopole" rest mass of the 
particle [6]. Using equations (4. 11), (4. 12) and (4.14) the equation (4.9) takes the following 
form 



1 DS"^ d / 

+ — + -r n — - A*"^ + ^A^° + — ^u„A''° 



2 ds ^ ds\ uO i V nO ^ -uO ^ 



(4.15) 



where round brackets in S^'^^u'*'^ denote symmetrization. Since M'^'^ is symmetric, the anti- 
symmetric part of the right-hand side of (4.9) vanishes. Hence one can immediately verify 

^ ^ - 4^ + ^ = A^^ - A^^ + ^(n-A^o - u'A'^^) (4.16) 
u" ds ds ds w' 

After some elementary calculations, this equation can be rewritten as 

-^+u''up^ nV— ^ = iA''^-A^'')+u''upiA^P-AP^)-u\p{A''P-AP'') (4.17) 

ds ds ds 

This equation should be compared to equation (5.3) in [3]. With the help of (4.15)-(4.17), 
we derive the following equations which are useful in simplifying (4.8) 

M^o = f mn^ + Up^^ - T^..^^"" - Up[A''P - A^^) j (4.18) 

r;.M- =r;, ( ,nu'u^ + + r;,rJ,s-„^ + r;,- (_^) ^^^^^^ 

By using equations (4.18) and (4.19), we can bring (4.8) into its final covariant form 



D / ^ DS^P 
— mu + u 



^ + ^SP'^u'^R"^,^ = £ (upiA'^P - AP")) + A'^ + y V^idx'^Snd^x 

(4.20) 



ds \ ds 

where R'^aufi curvature tensor. Also, note that it is straightforward to show 



n° J Va{dx"6nd^x ~ r^^(X)A"^ (4.21) 



Equation (4.20) should be compared to equation (5.7) in [3]. In general relativity, the right- 
hand side of (4.20) is zero [3]. 

It is important to note that the mass parameter m is not necessarily a constant of 
motion. To show this in more detail, let us define the "generalized momentum" p'^ as follows 

= mu" + ^—u^ + lupA^P""^ (4.22) 
ds 

where square brackets denote antisymmetrization. Using this equation and taking into ac- 
count that m = Uup'^ , it is straightforward to rewrite equations (4.17) and (4.20) as follows 

= 2 (pt'^n^l + aI'^^]) (4.23) 
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= As^^-u'^R\„^ + (a^ + r:;^A"/^) (4.24) 

where the dot denotes the covariant derivative with respect to the proper time, " ' " = 
D/ds. The main result of this section is embodied in the equations (4.23) and (4.24). These 
equations should be compared to the well-known Mathisson-Papapetrou equations for pole- 
dipole test particles in general relativity. It is obvious that if we assume that there is no 
coupling between matter and the gravitational fields, then A" and A"^ are zero and equations 
(4.23) and (4.24) recover the Mathisson-Papapetrou equations. 
Using equation (4.24) one can easily show 



m 



{js ^^'^"''''^ + 2^ZpA[H^ + (a- + r^^A"'^) (4.25) 



We can also define another mass parameter by = PaV"^ which in general is diff'erent from 
m. The mass parameter M sometimes called the "dynamical", "total" or "effective" rest 
mass of the test particle [6]. If we contract (4.23) with PuPx-, we get 



^ = -1?^- f^'^V) + - (a'^ + (4.26) 



Mm ds ( ~^ m 

Thus Ad is not also necessarily a constant of motion. 

It is also necessary to mention that just as in general relativity, equations (4.23) and 
(4.24) are insufficient to describe the motion of the pole-dipole particle. In fact, we have 14 
unknown variables (a;°, m° and S"^) and only 11 equations i.e. (4.24), (4.23) and u°'Ua = 1- 
Note that u°'Ua = 1 is the mere constraint due to the choice of parametrization of the 
trajectories. Thus, we need three more equations to make the model self-consistent. 

Let us recall the origin of this apparent inconsistency. The above procedure for deter- 
mining the motion of test particles consists of choosing a representative point (or a repre- 
sentative world line X^(t)) in the particle and taking the moments of equation (2.1) about 
that point. However, we have not specified how to choose this representative point in the 
test particle. In the other words, there is an arbitrary step in the above method, and it is 
the choosing of a representative point in the body. Thus, one can expect that if such a point 
is not uniquely specified, the motion of the representative point will not be fully determined. 
In general relativity, the proper specification of such a representative world line yields to a 
supplementary condition to Mathisson-Papapetrou equations which determines the particle's 
motion [22]. The two most widely used supplementary conditions are the Frenkel condition 
[23] (for massless particles) 

= (4.27) 
and the Tulczyjew condition [24] (for massive particles) 

5"% = (4.28) 

In general relativity, with the help of equations (4.25) and (4.26) one can immediately 
verify that m is constant when the Frenkel condition (4.27) is assumed, and the mass pa- 
rameter M is constant for the Tulczyjew condition (4.28). Here, we do not want to derive 
the corresponding supplementary equations which may yield to constant parameters m and 
M in the context of modified gravity theories. Because this may need the exact form of the 
nonminimal coupling term df^. However, we shall go into this issue for MOG in the next 
sections. 
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5 Test particle motion in MOG 

In this section, we consider the test particle motion in MOG, and we make some comments 
about this theory. MOG is an alternative theory of dark matter, and it is claimed in the liter- 
ature that this theory can address the dark matter problem in the spiral galaxies and clusters 
of galaxies [20, 25, 26]. MOG is a scalar-tensor-vector gravity theory which postulates, in 
addition to the metric tensor, three dynamical scalar gravitational fields G, fi and oj and also 
a dynamical massive four- vector gravitational field (j)^, see Appendix A for more detail. The 
vector field 4>^ is coupled universally to matter. As we mentioned before, such a coupling 
between matter and the vector field yields to a nonzero covariant divergence of the ordinary 
matter energy-momentum tensor, and consequently this coupling leads to a modification of 
the law of gravitation. Here by the law of gravitation, we mean the gravitational force law 
between two point particles in the weak-field approximation of the given modified gravity 
theory. We know that general relativity in the Newtonian limit coincides with Newtonian 
gravity, but this is not the case for other gravity theories. It should be emphasized that, 
for a modification to the law of gravitation, it is not necessary to construct a theory with 
nonminimal coupling terms between matter and gravitational fields. In fact, it is well-known 
that modified gravity theories which yield to a conserved energy-momentum tensor, can also 
provide a modification to the law of gravitation. For example see metric f{R) theory [27, 28]. 

In order to investigate the test particle equations of motion we need the exact form of 
the conservation law (2.1). We have found the conservation law in Appendix A. The result 
is (equation (A. 31)) 



where B^i, = V ^4>u — Vu(pfi, J"^ is a "fifth force" matter current defined by (A. 20) and the 
self-interaction potential W{(j)) is given by (A. 30). As we discussed, the right-hand side of 
this equation provides the modification of the law of gravitation. It seems that this theory 
can provide a wide range of modifications to the law of gravitation. In the other words, by 
choosing different self-interaction potential W{(j)) one can introduce different modifications. 
In [20], a special modification has been postulated (see equation (32) in [20]). This postulation 
can be considered as a special choice for W{(j)). However, we show that there is no potential 
W{(j)) corresponding to the extra force postulated in [20]. In the other words, we show 
that the test particle equation of motion postulated in [20] is not compatible with the field 
equations of MOG. Albeit, as we shall discuss, this issue may not change the main results of 



In this theory, the covariant derivatives of the vector field (j)^ do not appear in W{(j)). 
Thus, the most general form for W{<j)) can be given as 



where Cn are constant coefficients. For such a potential, it is straightforward to show that 
the second term on the right-hand side of (5.1) vanishes. Finally, the conservation law of the 
ordinary matter energy- momentum tensor in MOG reads 




(5.1) 



MOG. 




(5.2) 



n 



(5.3) 
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Comparing this equation to (2.1), we get Sf^ = ^—gB^SJ'^. In the presence of the test 
particle, the matter current density is given by J'^ = + 5J^ . On the other hand, since 
the particle is moving outside of the massive bodies is zero near and inside of the test 
particle. Furthermore, since it is assumed that the matter current density is conserved (i.e. 
Vq-J"^ = 0), one can immediately infer that VctdJ" = 0. Consequently it is straightforward 
to show that Q = J yJ—g5J^d^x is a time independent quantity. Hence, we may introduce 
the matter current density 5 J" as follows (see [29]) 

Sr = -%= 1 6^{x- X)dx" (5.4) 
J 

Therefore, it is easy to show that A'^ is given by 

A'^ = gs^ (5.5) 

where = —Q is the test particle's fifth force charge. In this theory, the fifth force charge 
is assumed to be proportional to the particle's mass m, such that = nm. Finally, using 
equation (3.17) the equation of motion of the single-pole particle is given by 

mu^ = q5B^u°' (5.6) 

Taking into account equation (3.14) it is clear that as in general relativity m is a constant 
of motion, and can be considered as the rest mass of the test particle. Equation of motion 
(5.6) is different from what postulated in [20] (equation (31)). In fact as it is obvious from 
(5.6) the scalar field uj{x), and its derivatives do not appear in the equation of motion. As 
we mentioned before the vector field (j)'^ is coupled to ordinary matter, thus we expect that 
the vector field or its derivatives exist in the equation of motion. On the other hand, there is 
a coupling between (/>^ and u}{x), see equation (A. 3), and consequently one may expect that 
the scalar field u!{x) should appear in the equation of motion. However, as it is clear from 
(5.6), it is not the case in this theory. 

In [20], it has been assumed that, for the vacuum solution around a spherically sym- 
metric point-like mass, uj{x) is nearly constant. This assumption is also consistent with the 
numerical solution of the vacuum field equations [30]. With this unnecessary assumption, 
the postulated equation of motion coincides with (5.6), and the main result of this theory 
related to the dark matter problem will not change. 

Now let us consider the spinning particles equations of motion in MOG. It is straight- 
forward to verify that A'^^ is zero in MOG. Therefore, the equations of motion (4.23) and 
(4.24) can be simplified as follows 

S"^ = Ip^^u^^ (5.7) 

= _ is^^n'^i?-^,^ + q.B-'^u'' (5.8) 

where p'^ = mu'^ + S'^'^Ua. Equation (5.7) is the same as in general relativity, but there exists 
the extra term q^B'^^u" on the right-hand side of (5.8) relative to Mathisson-Papapetrou 
equations. It is obvious from equations (5.7) and (5.8), which describe the motion of electri- 
cally neutral pole-dipole particles in MOG, that there is a similarity between these equations 
and the Dixon-Souriau equations which determine the motion of charged pole-dipole particles 
in the presence of electromagnetic fields in the context of general relativity. Mathematically, 
the origin of this similarity lies in equation (5.3). In fact, using Einstein-Maxwell equations 
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one can easily show that the covariant divergence of the ordinary matter energy-momentum 
tensor is given by 

V^T^'^ = Fjr (5.9) 

where F^jy is the field strength tensor, = {p, , , J'^) is the current four- vector and p 
is the electric charge density. The close similarity between equations (5.3) and (5.9) yields 
to a similarity between equations of motion of test particles in MOG and Dixon-Souriau 
equations. However, one should note that equations (5.7) and (5.8) are not exactly similar 
to the Dixon-Souriau equations. In fact, mathematically, they are similar to the equations 
of motion of a charged spinning particle without magnetic moment, see equation (9) in [13] . 



6 Conserved quantities 

It is well known that the symmetries of the background space-time may yield to the existence 
of conserved quantities. For example, in general relativity, for a charged test particle moving 
under the influence of electromagnetic field, if the electromagnetic field satisfies some consis- 
tency conditions, then the space-time symmetries guarantee the existence of some conserved 
quantities. In this section, we consider the conserved quantities related to the existence of 
the Killing vectors and the Killing tensors. We recall that a Killing vector and a totally 
symmetric rank-n Killing tensor Kp^ j^^ satisfy 

V(„C/3) = 

and 

V(„i^^,..^„) = (6.2) 

Let us first consider the conserved quantities for the motion of a single-pole particle. The 
equation of motion is given by (5.6). If we assume that the background geometry possesses 
a Killing vector then one can straightforwardly show 

^ (CaP° + (Zs^a^ ) = q^u^'hckt. (6.3) 
where L^(j)^ is the Lie derivative of the vector field with respect to ^. Thus the quantity 

Qi = iaP'' + q^io^r (6.4) 

is conserved if the vector field (j)^ is Lie-conserved along i.e. L^^cf)^ = 0. Similarly, one can 
verify that the component of the velocity along ^ 

= (6.5) 

is conserved if 

iaB% = (6.6) 

This condition can be generalized to quantities nonlinear in velocities if the space-time has 
appropriate Killing tensors. Let -ftT/^j.../^^ be a rank-n Killing tensor satisfying equation (6.2). 
Then it is straightforward to show 
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Thus, Qk is conserved if 

K\p,...p^_Bp^y = (6.8) 

Since the metric g^p is a rank two Kilhng tensor satisfying triviahy the condition (6.8), we 
find conservation of the norm of the test particle's velocity i.e. g^^yU^u'^. 

Now let us consider the conserved quantities for the motion of a spinning test particle 
in MOG. Using equations (4.25), (4.26) and (5.7), we find 



D 

ds 



m = — iS^'^up'^ ila 

where S is the magnitude of the spin defined by = ^Sq^/S"^. Thus, m is a constant of the 
motion when the Frenkel condition (4.27) is assumed and accordingly M is constant when the 
Tulczyjew condition (4.28) is assumed. Also, the magnitude S of the spin is a constant of the 
motion for both supplementary conditions. Furthermore, if we assume that the background 
geometry possesses a killing vector then using equations (5.7) and (5.8) we get 



D 

ds 



(CaP" + g5Car + ^S^'^V^C/?) = q5U^'L^^^. (6.10) 



Thus for a spinning test particle the quantity 

Qi = ic.p'' + q^iar + \s^^'^o.ip (6.11) 

is conserved if the vector field (f)^ is Lie-conserved along ^. 

A Covariant divergence of T^" in MOG 

The action of the theory is given by 

S = Scrav + Sif, + + + + Sm (A.l) 
where 5m is the action for ordinary matter and 



^d^xio(^^B^^''B^, + V44>,g^,,^l^)^ (A.3) 

S^ = JV^ d^x Q^^'^V^xV.x - V^U?j (A.4) 

S^ = jV^ d^x fix) Q^'^'v^vv.v - Ka(V')) (A.5) 

S^ = d^x fix) Q^^'^V^L.V.w - K;(o.)) (A.6) 
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where R is the Ricci scalar, Bf^^, = V — Vycf)^ and the scalar fields x ^-^id ip are related to 
the scalar fields G and fi of the original paper [20] by 

/(X) = Y = ^, V' = ln// (A.7) 
also and are related to Vq and V^j by 

Vr V„ 

yx = f, y^ = f, (A.8) 

In this theory, it has been assumed that gy^^ = 0. In other words, the covariant derivatives 
of do not appear in V^. 

The total energy-momentum tensor is given by 

rptotal _ rp _i_rp4) , rpX _|_ T^V' _|_ I A Ol\ 

where 

— 2 6Sm rnO ~^ IK iri\ 

and Q can be cp, Xj V' or u. Using equations (A.2)-(A.6) and (A. 10), one can easily verify 

= ^ ( - 9,u [\b'^''B,„ + V^^ + 2^) (A.ll) 

T;^. = - (v^xV.x - (^V,xV°x - ^x) ) (A- 12) 

= -fix) (v,,VV.V - Qv«^V> - K/,)) (A.13) 

= -fix) (v^^V.w - Qv«6.V"^ - (A.14) 

Furthermore, variation of (A.l) with respect to g^'^ , (pai Xi "0 'ind u yields to the following 
field equations respectively 

fix)R^u - lfix)Rgt.u = (v^v. - 5^.n)/(x) + SttT^"*'^' (a.15) 

wV^5°^ + V^a;B"^ + ^^-/ = -J" (A.16) 

ax + V; = ^ + fix) ^^^'^(V^^V.V' + V^^V^o;) - (V^ + K.)^ (A.17) 

f{x)iD^ + V;) = -/'(x)V,xV^V' - (A.18) 

/(X)(nu; + V:,) = -/'(x)V^xV^^ - ^5^"i3p. - (A.19) 
where prime stands for H'{x) = ^ and J" is a "fifth force" matter current defined as 

J" = --L^ (A.20) 
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Taking the covariant divergence on both sides of equation (A. 15) yields 



V^r^. = "T^^V.x - (V%t + V%^, + V%t + V%"J (A.21) 



Note that on purely geometrical grounds, (DV^ — V^D)/ = Rf_iuV^ f and V^G^^ = 0, where 
/ is an arbitrary scalar function and G^i, is the Einstein tensor. On the other hand, taking the 
covariant divergence of equations (A.11)-(A.14) and using the field equations (A.16)-(A.19), 
we find 

V^T/, = J" - (^\bp^B,^ + V^^ + 2V'' (w^) - ^^uV^ + ^B,^^ (A.22) 

V%^, = -^^V-X - /'(x)V.x ^^"^(VaV'V^V + V^^V^a;) - (F^ + V^)j (A.23) 

v^r/f. = fix) Qv„^v„^ - K/,^ V,x + V.V' (A.24) 

V%", = fix) Qv,t^V«w - K.) V,x + i^^BP'^Bp^ + (A.25) 
In equation (A.22) we used the following equation 

B^'V^B^a - \v,{BP''Bp„) = (A.26) 

Equation (A.26) can be easily verified by using the definition of B^^^. Now substituting 
equations (A.22)-(A.25) into equation (A.21) we obtain 

V^T,, = B^^r + (c^V«</'.|^ - 2V^ (^^)) (A-27) 

On the other hand, taking the covariant divergence of equation (A. 16) and assuming Vq, J" = 
0, one gets 

V„ fw^l = (A.28) 



With the help of this equation we can rewrite (A. 27) as 

V^T,, = B^^r + (5Ma^<A.|^ - 2a;p;) (A.29) 
Finally, we choose the following form for the potential (see equation (20) in [20]) 

V^ = -\e^'^(p''(Pa + W{(p) (A.30) 

where W{(j)) is the vector field cj)"' self-interaction contribution. Substituting this potential 
into (A.30) we get 
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